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The idea is considered that a classical non-Abelian gauge field can be considered as a dark matter 
candidate. It is shown that Yang-Mills equations have solutions with such distribution of the mass 
density that allows to describe a rotational curve of spiral galaxies. The conditions necessary for 
such consideration are considered. One parameter is estimated from Yang-Mills equations and from 
astrophysical observations (stars rotation curve). The agreement is to within 1%. 
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I. INTRODUCTION 

The electromagnetic field behaves at different conditions either as pure quantum or almost as a classical field. But 
for non-Abelian gauge fields we know only quantum manifestations of these fields. Here we would like to present the 
idea that classical non-Abelian gauge fields can become visible as DM. 

The problem of the DM nature is now one of the most fundamental problems in modern physics. The reason is 
to explain the invisibility of DM it is necessary to involve such kinds of an exotic matter as, for example, WIMP, 
supersymmctric particles etc. 

In this article we offer a new model of DM. In the studied model the DM is an interior of an abortive singularity 
or black hole. We understand this suggestion as follows. 

In Ref 's 01 [1| it is shown that the energy density of spherically symmetric solutions of Yang-Mills equations ma' 
have non-standard behavior (behind an exception, of course, 't Hooft - Polyakov monopole). For example, in Ref. 
it is shown that the energy density is weakly decreasing. It is natural that such distribution of the matter will give 
diverging solutions by switching on the gravity. For some boundary conditions at the center it will be an interior of 
a non-Abelian back hole and for other values of the boundary conditions a singularity located at some distance from 
the center will appear. But one very interesting subtlety exists here. This solutions arc classical ones. The careful 
analysis shows that at the infinity these solutions are strongly oscillating in the space (gravity is switched off) and 
the period of the oscillations increases by moving away from the centre. It is obvious that on some distance from 
the center the quantum fluctuations of gauge field on the period distance of oscillations become comparable with a 
field magnitude. It means that at this distance the gauge field becomes quantum one. The problem in calculating of 
such distribution of the non-Abelian gauge field is that the gauge field becomes essentially non-pcrturbativc and it is 
impossible to apply the perturbative Fcynman diagram techniques to its description. 
^ . The idea presented here is that the radius where the gauge field becomes quantum can be less than the radius of a 
' singular or an event horizon corresponding to a non-Abelian black hole. In this case this object docs become neither 
, a singularity nor a black hole (abortive singularity/black hole). It looks as follows: in space there is a sphere filled 
' with a classical gauge field and on border of this sphere the non-Abelian gauge field becomes quantum one and the 
mass contained in this sphere is not enough for formation of a singularity or an event horizon. The color gauge field 
in the sphere does not interact with the elementary particles because the particles are colorless and consequently this 
classical field is invisible and can be applied as a candidate of DM. 

In fact in this paper we consider the idea that in a non-Abelian gauge theory may exist such space distribution 
of a gauge field that classical and quantum phases exist simultaneously but spatially separated. In this case the 
jump condition from classical phase to quantum one is the strong oscillations of classical non-Abelian field. These 
oscillations leads to the fact that quantum fluctuations become essential at a distance of oscillation period. In other 
words the space is filled with a non-perturbative gauge vacuum where there exist defects filled with the classical gauge 
field. The galaxies are located in these defects and the classical gauge field is the DM. 
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II. THE INTERIOR OF AN EINSTEIN- YANG-MILLS SINGULARITY/BLACK HOLE 

In this section we would like to show that usually the SU(3) gauge field distribution leads to a singular spacetime. 
Wc use the following metric 
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where t, r, 6, ip are usual spherical coordinates. Substituting metric ([T]) and SU(3) gauge potential (jAip - (|A6P into 
Einstein- Yang-Mills equations 
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and Yang-Mills equations 
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where x is the gravitational constant, a ~ 1, 2, • • • , 8 is the color index and xi''') — h{r) — 0. The analytical solution 
does not exist and we search the numerical solutions for this equations set. 

We search for the solution inside of the non-Abelian singularity/black hole. It means that wc should start the 
solution from the point r = 0. The analytical solution close to the center is 
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The numerical solution can not start from the point r = since Eq's ©-([T]) have terms l/r^. Consequently we should 
start from the point x = S <^ rj^fyt. The boundary conditions are 
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We present the profiles of the functions g""^ = (1 
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),i^{r), v{r),w{r) in Fig's[T]-[3] We see that at a point r = vh 
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the function g''^ (r^) = 1 — = 



0. The energy density is 
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and its profile in Fig. [4] is presented. Consequently the numerical analysis have shown us that at the point rn where 
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FIG. 3: The profile of the function ■w{r). 



FIG. 4: The profile of the function energy density. 



1 — M{rH)/f'H — there is a real or coordinate singularity (for some boundary conditions it will be an event horizon). 
From Fig. [2] we see that the function v{r) is oscillating function and the period of oscillations increases with the 
moving away from the center. At some distance from the center the quantum fluctuations on the distance of period 
will be comparable with the field magnitude. In this case the classical gauge field becomes quantum one. If it occurs 
on the distance < th then we will have a ball filled with the classical gauge field inside quantum gauge field. In 
Ref. [l3| is claimed that the non-perturbative quantized SU(3) gauge field can be described using two scalar fields. 
These fields decrease very quickly (exponentially) to a ground state. If this so then the situation looks as follows: in 
the space there is a ball filled with the classical gauge field with weakly decreasing energy density. At some distance 
Tg from the center the field becomes quantum one and the energy density very quickly (exponentially) decreases to 
a ground state (non-perturbative vacuum). The idea presented here is that the ball can be considered as the DM as 
usual elementary particles are colorless and consequently can not interact with the color classical SU(3) gauge fields. 

Let us note that non-Abelian black holes and particlelike solutions exist only with a special choice of the boundary 
conditions on the event horizon or at the center. We have to underline that in contrast with the non-Abelian black 
holes and particlelike solutions we consider the solutions with arbitrary boundary conditions. 

In the following sections we will estimate the radius ru of a real/coordinate singularity and the radius Vq where the 
transition from the classical phase to quantum one occurs. 
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III. THE ESTIMATION OF THE RADIUS OF REAL OR COORDINATE SINGULARITY 



Unfortunately the solution presented in previous section is numerical one that does not allow us to calculate the 
sphere where the above mentioned singularity (real or coordinate) is. We will estimate the radius of a singular sphere 
(using Newton gravity) in the following way. The radius is estimated as the place where the Newtonian gravitational 
potential becomes so strong that a test particle should have a velocity of light to get from this point on infinity. We 
test this method for the Schwazschild black hole solution. The energy conservation law tells us 
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here G is the Newtonian constant; the first and second terms are kinetic and potential energies of a test particle with 
the mass m and M is the mass of a singularity or a black hole. As a result we have the radius of a singularity or an 
event horizon 
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that absolutely precisely coincides with the event horizon radius calculated in general relativity. 

Now we will try to estimate the radius of the singularity or the event horizon for the classical distribution of the 
SU(3) gauge fields (|A8P - (IAlip . We will work in Newton gravity because the values V2 and ws are small enough, i.e. 
the magnitudes of the SU(3) gauge fields close to the center are small enough. In this case Yang-Mills equations in a 
flat space are 
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here the dimensionless radius x = r/ro is introduced and rg is an arbitrary constant. In fact Eq's (|20[) (|2ip arc Eq's 
([5]) (O in Minkowski spacetime. The asymptotical behavior x- ^ 1 of the solution is Q 
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where v' = dv/dx,w' = dw/dx. Asymptotically the dimensionless energy density is 
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The numerical analysis shows that the functions v{r) and w{r) quickly attain an asymptotic form. Therefore we will 
use (|26p to estimate the mass (r) under radius r 
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that leads to 
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where = ^^7^ is the dimensionless couphng constant in the SU(3) gauge theory; 1/g is the analog of the electric 



charge in electrodynamics; Ipi = 



10 m is the Planck length. The numerical factor 
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and consequently 
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Using the parameters rg and a which are calculated in Appendixes |B] and [C] we have 



rH « lO^^m 



(32) 



that is much more that the galaxy radius. 

Now we can define the natural choice of the parameter rp. At the origin Eq's (j20p (j21[) have the solution that can 
be presented as series (|11|) 
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Therefore the natural choice of the parameter rg is 
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IV. THE TRANSITION FROM CLASSICAL PHASE TO QUANTUM ONE 

The idea considered above mentioned solution of the classical Yang-Mills equations can not be extended up to 
infinity because in some place the space oscillations of the classical gauge field becomes so strong that non-perturbative 
quantum effects should be taken into account. In this section we try to estimate the radius where it can happen. 

Following on the Heisenberg uncertainty principle 
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here AF/^ is a quantum fluctuation of color electric field F/J ; A is a quantum fluctuation of color electric potential 
A°*; Ay is the volume where the quantum fluctuations AF^^ and Aj4'" takes place; a = 1, 2 • • • 8 is the color index; 
i — 1,2,3 is the space index. 
For the ansatz (|A8l) - (|Alip 
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We introduce the physical component of the F^g 



2 . vw 
— smO—^. 

g 



(38) 



6 



To an accuracy of a numerical factor the fluctuations of the SU(3) color electric field are 

AF^eK-^{Avw + v Aw). (39) 

For the ansatz ((X8|) - (|Alip 
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Introducing the physical components of the gauge potential Ag-^'^-'^'^ 
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The volume AV is 

Ay = 47rr2Ar. (44) 
The period of space oscillations by r ^ rg can be defined in the following way 
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We suppose that the place where the SU(3) classical color field becomes quantum one is defined as the place where 
the quantum fluctuations in the volume AV = Anr^Ar with 
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of the corresponding field becomes comparable with magnitude of these fields 

Av « V, Aw ~ w (47) 
Substituting of Eq's (gS]), ^ , ^ and ^ into Eq. ^ we obtain 
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where — -^j- is the dimensionless coupling constant that is similar to the fine structure constant in quantum 

2 2 

electrodynamics a = j- . In quantum chromodynamics /3 = 1/g' > 1. If we choose 1/g' « 1 and from Fig. [2] we take 
A « 0.4 we see that 

2 

« 6.25 (49) 



that is comparable with 2tt w 6.28. 

Thus in this section we have shown that if the condition (j48p is true then at some distance from the center the 
transition from the classical phase to quantum one occurs. Unfortunately the rough estimation presented in this 
section docs not allow us to calculate the radius where such transition takes place. For the exact evaluation of the 
place where such transition happens it is necessary to have non-perturbative quantization methods which arc missing 
at the moment. 



V. THE ROTATION CURVE OF YANG-MILLS COLORED DARK MATTER 



In this section we would like to show that the solution of Eq's (I^D)) (PT|) really has such mass density distribution 
that it is in a good agreement with the Universal Curve Rotation. 
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A. Numerical investigation 



111 this subsection we present the typical numerical solution of Eq's ((20)l pT|) . For the numerical investigation 
we have to start from the point x ^ S <^ 1. Here we have approximate solution (|33p (j34p and now we choose the 
parameter vq as 
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The typical behavior of functions v{x) and wlx) is presented in Fig.[5l 
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FIG. 5: The profile of functions w{x),v{x),e{x), V2 = 
-0.1, W3 = 1. 




FIG. 6: The profile of the dimensionless rotation curve 



The mass density p{r) is 
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where /^(a;) = £{x) and e(a;) is given in Eq. (j25p . The profile of the dimensionless energy density e{x) in Fig. [5] is 
presented. 

The rotation curve is defined as 
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where A4{x) is the dimensionless mass of the color fields A° inside the sphere of radius r = xrg, g' = g'^ch/Aw is 
the dimensionless coupling constant, G is the Newton gravitational constant. The parameter a w 1.31 can be found 
using fitting of functions w{x) or e{x), for details see Appendix IbI 



B. The comparison with a Universal Rotation Curve of spiral galaxies 



In Ref. [ill a Universal Rotation Curve of spiral galaxies is offered that describes any rotation curve at any radius 
with a very small cosmic variance 
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where Ropt = 3.2 Rjj is the optical radius and Rd is the disc exponential length-scale; X = r/ Ropt] L is the luminosity. 
We would like to compare the rotation curve for the color fields ([52]) with the Universal Rotation Curve ((53|) where, 
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for example, L/L^, = 1 
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For the DM the Universal Rotation Curve is 
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The profiles of y(7i?c(^), V^j^^{X),Vl,^j{X) in Fig. [7] are presented (Vg^^ is the rotation curve for the DM, V^m(^) 
is the rotation curve for the light matter). 
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FIG. 7: The profiles of dimensionless rotation curves for 
the light and dark matter [lH . 



FIG. 8: The comparison of DM rotation curve (|55p 
(curve 1) with the rotation curve (|6ip (curve 2) for the 
SU(3) classical color field a f» 1.31, 5' = 

1, Ropt = 2OKP3 = 6 * 10"m, Vopt = WOKm/s. 



At the center r w the approximate solution has the form p3|) (p4l) and the mass density (j5ip approximately is 
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Consequently the rotation curve will be 
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Far away from the center the dimensionless energy density ecx^{x) is presented in Eq. ((26)) and in this case we can 
estimate the values of square of speed in the following way 
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The numerical value of is defined near to the center of galaxy where according to Eq. ([60)) the difference £Qo{x)—e{x) 
is maximal. Thus the asymptotical behavior of the rotation curve for the domain filled with the SU(3) gauge field is 

- Vl (61) 

In Fig. [5] the profiles of the Universal Rotational Curve ([55)1 and fitting curve ([ST)) are presented. The value of 
parameter a is given from section IV Al The details of fitting 

ro w 2.01 • 10-1^ cm, (62) 
Vb w 32.25 Km • s'^ (63) 

in Appendix IC) are presented. One can see that the biggest disagreement is near to the center since close to the center 
the fitting curve have to be ([57)) not ([6T)) . 
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C. Transition to non-perturbative quantized phase 

The energy density ()25p gives us an infinite total mass. How can we avoid this problem ? In Section [TV] we brought 
forward arguments that the gauge field inside of some region is in classical phase and outside the region is in 
quantum phase. We think that it is the manifestation of the fact that the gauge field is the strongly interacting field 
and the quantization of this field should be carried out using a non-perturbative technique. In fact in this paper we 
would like to show that non-perturbative quantized fields can be spatially distributed in such a way that classical and 
quantum phases exist simultaneously. 

Unfortunately up to now we do not have any exact non-perturbative technique for the quantization. In this section 
we want to describe briefly approximate non-perturbative technique based on the Heisenberg approach for non- 
perturbative quantization of a nonlinear spinor field (for details, see [loj). 

In section HV] we have shown that at some distance from the center the classical phase changes on non-perturbative 
quantum phase. In Ref. [l3| it is shown that two scalar fields may describe a non-perturbative quantized gauge field. 
Briefly it can be shown by the following way. In quantizing strongly interacting SU(3) gauge fields - via Heisenberg's 
non-perturbative method [l3| one first replaces the classical fields by field operators A^. This yields the 

following differential equations for the operators 

d^^T^'"' = 0. (64) 

where fj,,iy — 0,1,2,3; B = 1, 2, • • • , 8 are SU(3) color indices. These nonlinear equations for the field operators of 
the nonlinear quantum fields can be used to determine expectation values for the field operators A^. One problem in 
using these equations in order to obtain expectation values like {A^), is that these equations involve not only powers 
or derivatives of {A^) {i.e. terms like da{A^) or dadp{A^)) but also contain terms like Q^^ = {A^A"^). Starting 
with Eq. ([64p one can generate an operator differential equation for the product A^A'^' consequently allowing the 
determination of the Green's function Gj^^ 

q) = 0. (65) 

However this equation will in it's turn contain other, higher order Green's functions. Repeating these steps leads to 
an infinite set of equations connecting Green's functions of ever increasing order. This construction, leading to an 
infinite set of coupled, differential equations, does not have an exact, analytical solution and so must be handled using 
some approximation. The basic approach in this case is to give some physically reasonable scheme for cutting off the 
infinite set of equations for the Green's functions. Using some assumptions and approximations on 2- and 4-points 
Green's functions one can reduce the initial SU(3) Lagrangian to an effective Lagrangian describing two interacting 
scalar fields (for details see Ref. [l^l)- The scalar fields and x which are under discussion here appear in the 
following way. We assume that in the first approximation two points Green's functions can be calculated as follows 

{A'^(x)A]{y)) = -rfc/"P"7''^"x'"(x)x"(y), (66) 
(Ag(x)A^(y)) « {A'i{x)A]{y)) (67) 
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where € SU{2) C SU{3), a = 1, 2, 3; m = 4, 5, 6, 7, 8; i, j = 1, 2, 3 are spatial indices. And 
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where A™ e SU{3)/ SU{2). The 4-points Green's functions are a bihnear combination of 2-points Green's functions 
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here Ai,2,Mi,2 are some constants. The assumptions (|66l) - (f7T|) allows us to average the SU(3) Lagrangian 
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and bring it to the form 
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with the field equations 
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Let us consider the spherically symmetric case — k(f){r),x"^ = kx(r) where k is some constant. In this case the 
field equations are 
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It is easy to see that asymptotically the solution has the form 
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where 0cx3 , Xoo , ^'ij are constants. We think that this solution describes the non-perturbative quantized SU(3) gauge 
field after the transition from classical phase to quantum one occurs. 

The main point of this consideration is that the non-perturbative quantized gauge field decreases very quickly 
(exponentially) after transition to the quantum phase and consequently the total mass becomes finite one. 
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VI. INVISIBILITY OF COLORED DARK MATTER 

For the detection of DM (in the context of the model presented here) we can use only colored particles which can 
interact with color gauge fields. The equations describing the motion of a colored particle are Wong's equations 



m 



where x^{s) is the 4D trajectory of the particle with the mass m. is the color components of color charge of the 

particle, (Q'^) = ~ const. The ordinary elementary particles are colorless and consequently do not interact with 
the color DM. 

Only 't-Hooft - Polyakov monopoles and dyons may interact with colored DM and can be used for the detection of 
the colored DM. Another possibility of the interaction of an elementary particle with colored DM is the interaction 
between external color field (DM) and a color electric and/or magnetic dipole or quadrupole of this elementary particle. 



VII. CONCLUSIONS 



In this paper we have considered the idea that the problem of DM probably can be connected with the problem 
of non-perturbative quantization of strongly interacting fields. It allows us to connect one problem in macroscopical 
physics (the problem of DM) with another problem in microscopical physics (confinement problem on quantum 
chromodynamics). In this connection R. Kolb in Ref. [l3| write: "Dark matter and dark energy are two of the binding 
cords I will use to illustrate how collaborations of astronomers and high energy physicists on large astronomical 
projects can be good for astronomy, and how discoveries in astronomy can guide high-energy physicists in their quest 
for understanding nature on the smallest scales." . 

The features of the model of colored DM presented here are: 

• The estimation of the parameter a from the Yang-Mills equations (see Appendix [B]) and from the astrophysical 
(see Appendix [C]) point of view is in agreement with remarkable accuracy. 

• On the background of a non-perturbative vacuum of non-Abclian gauge field there exists a bubble of the same 
classical gauge field. The classical non-Abelian gauge field filled the bubble is colored DM. 

• Spherically symmetric classical solutions of the Yang-Mills equations have weakly decreasing mass density dis- 
tribution leading to a good agreement with the Universal Rotation Curve. 

• The distribution of these classical gauge field is that on some distance from the center the transition from 
classical phase to quantum occurs. 

The problems for the future investigations are: 

• The theoretical estimation of Vq parameter and comparing it with the fitted value (j63p . 

• The fitting of the Universal Rotational Curve ([55]) using the function joining ([ST)) and (|61|) . 

• The most important problem in the model presented here is the calculation of gauge field distribution using a 
non-perturbative quantization technique. 

• The search for possibility of the classical gauge fields detection . 
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APPENDIX A: ANSATZ FOR SU(3) GAUGE POTENTIAL IN MINKOWSKI SPACETIME 

We consider the classical SUTS) Yang-Mills gauge field and use the following ansatz for the SU{2) e SUii) 
components of the gauge field 

Al = -2^xW, Al = 2^xW, Al ^ -2^xW, (Al) 

gj.^ gj,Z gj.^ 

A\ = 2^[Mr) + l], (A2) 

g^l 

A\ - -2^ [%) + !], (A3) 
= 2^ \h{r) + 1] (A4) 

where ^)f^, = df,A^ - di,A^ + gf^^'^A^A^ is the field strength tensor; f^^^ are the SU(3) structural constants; 
A, B,C = 1,2, ■■ ■ ,8 are color indices; g is the coupling constant and 

{A,)^p = 2{e,,^x^' + e,spx") ^v{r), (A6) 

for the coset components belonging to the coset space SU{3)/ SU(2); z = 1, 2, 3 are space indices; e^j^ is the absolutely 
antisymmetric Levi-Civita tensor; the functions x{i')Th,{i')Tw{''^)T''^{i') ^-rc unknown functions. The coset components 
iA^)^p in the matrix form are written as 

a=l,3,4,6,8 

where = ^ are the SU(3) generators, are the Gell-Mann matrices. 
This ansatz in the spherical coordinate system is 

Af — {it;(r) sin^ 6'sin(2iy9); — 20(r)cos0; w(r) sin^ 6* cos(2iy9); w(r) sin(20) cos 0; 

20(r)sin6'sin^; u;(r) sin(26') sin </>; - (f>{r) sm{20); w{r) ^ +3cos(26>) | ^ ^^^^ 

2v3 J 

A^ = 0; (A9) 
Ag = {-2w(r) cos (2(^) sin0; 0; 2t;(r) sin sin (2(^) ; 2t;(r) cos sin v?; 2 [1 + h{r)] cos (p; 

-2w(r)cos6icos(^; 2 [1 + ft.(r)] sin </?; 0}; (AlO) 
= {!'(?') sin(26')sin(2(p); -2[1 + h{r)] sin 9; v{r) sin{29) cos{2ip); 2v (r) cos{20) cos ip; 

— 2[1 + /i(r)] cos6'sincp; 2z;(r) cos(26') sin (^; 2[1 + ft.(r)] cos 61 cos (^; V3w(r)sin(20)|. (All) 



APPENDIX B: FITTING OF PARAMETER a 



For the estimation of the parameter a we use the functions w{x) and e{x) given from the numerical solution of 
Eq's (PO)) (|2ip (see also Fig. O in the region 20 < a; < 100. The fitting functions are the asymptotical form of the 
function w{x) (sec Eq. ((23| ') 

'w{x) w ax" + Wo (Bl) 
and for the control we use the asymptotical form of the function e{x) (see Eq. ([26| ) 

£oo(a;) « ^a^ (a - 1) (3a ~ 1) x""-^ + Eq (B2) 

here Wq and Eq are systematical errors in the consequence of ignoring of other components in the asymptotical 
decompositions (|B1[) and (|B2p . 



13 



The fitting carried out using MATHEMATICA package. The fitting parameters are: a and either wq or Eq. The 
result of fitting is 

a « 1.31077, Wo « 2.15369 (B3) 

for fitting w{x) and 

a « 1.31995, So « -0.0012265 (B4) 

for fitting e{x). 

APPENDIX C: FITTING OF ROTATIONAL CURVE OF GAUGE FIELD 

For the fitting of the rotational curve ((6T|) we use the data from the Universal Rotational Curve ([55|) . The fitting 
equation is equation (|6ip in the form 



XL 

^opt 



Ax^ + C, (CI) 



^ ^ ^ ^^(a-l)(3a-l) 1^ /IpiV rRopt_Y'~' (C2) 



2 a^(a- l)(3a- 1) 1 / lpi_y ( Ropt 

3 2a - 1 V^p, \ro) \ ro 

B = 2a -2, (C3) 

C = (C4) 



where the fitted parameters are A, B, C. The fitting is carried out using MATHEMATICA package. The fitting 

a = 1.31954, ro w 2.01 • 10"^^ cm, K) ~ 32.25 Km • s~\ (C5) 



parameters are: ro and . The result of fitting is 
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